Research  and  Development  Technical  Report 
SLCET-TR-91-13 


A  Computer  Simulation  of  an  Adaptive 
Noise  Canceler  with  a  Single  Input 

Stuart  D.  Albert 

Electronics  Technology  and  Devices  Laboratory 


June  1991 


Approved  for  public  release. 
Distribution  is  unlimited. 


91-12714 

_ a-ai 


U.  S.  ARMY  LABORATORY  COMMAND 
Electronics  Technology  and  Devices  Laboratory 
Fort  Monmouth,  NJ  07703-5601 


NOTICES 


Disclaimers 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  desig¬ 
nated  by  other  authorized  documents. 

The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 
referenced  herein. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704  0188 


*JUK  ttoomn}  bufd*n  tot  urn  toflecnon  o*  m(o'm»t>on  it  nuimtrd  to  «v«r» g«  i  hour  o*  rhoont*.  including  th«  t,m*  tot  rrvwwtng  inwruroom.  Marching  tinting  d»u  tuu «*». 
gathering  and  maintaining  <h«  data  nttdtd.  tod  completing  ind  reviewing  iht  collection  of  information.  Send  comment,  regarding  thn  burden  etUmete  or  any  othef  eipect  of  thn 
collection  of  information.  including  wggetnont  tor  reducing  ihn  burden,  to  Wethmgton  Headquarter,  Servicet.  Directorate  for  information  operetions  end  Reports  tllS  letlerton 
Dews  Highway.  Suite  1104.  Arlington.  VA  22J03-TJ0J  end  to  the  Office  of  Management  end  Rodget.  Peoecwort  Reduction  Rro|ect  (0704-g  tM).  Washington  DC  JOSOJ 


1.  AGENCY  USE  ONLY  (Let ve  blink)  I  2.  REPORT  DATE 

(  June  1991 


4.  TITLE  AND  SUBTITLE 

A  COMPUTER  SIMULATION  OF  AN  ADAPTIVE  NOISE  CANCELER 
WITH  A  SINGLE  INPUT 


«.  AUTHOR(S) 

Stuart  D.  Albert 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSEES) 

US  Army  Laboratory  Command  (LABCOM) 

Electronics  Technology  and  Devices  Laboratory  (ETDL) 
ATTN:  SLCET-M 

Fort  Monmouth,  NJ  07703-5601 


9.  SPONSORING /MONITORING  AGENCY  NAME($)  AND  ADORESS(ES) 


J.  REPORT  TYPE  AND  DATES  COVERED 

Technical  Report:  1988-1991 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


SLCET-TR-91-13 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


1 J  ABSTRACT  (Minimum  200  word!) 

I 

j  A  description  of  an  adaptive  noise  canceler  using  Widrows1  LMS  algorithm  is  pre- 
j  sented.  A  computer  simulation  of  canceler  performance  (adaptive  convergence  time 
]  and  frequency  transfer  function)  was  written  (for  use  as  a  design  tool).  The 
simulations,  assumptions,  and  input  parameters  are  described  in  detail.  The  simu- 
;  lation  is  used  in  a  design  example  to  predict  the  performance  of  an  adaptive  noise 
|  canceler  in  the  simultaneous  presence  of  both  strong  and  weak  narrow-band  signals 
j  1  :•  cosited  frequency  hopping  radio  scenario). 

,  ■  7he  basis  cf  the  simulation  results,  it  is  concluded  that  the  simulation  is 

’ele  for  use  as  an  adaptive  noise  canceler  design  tool;  i.e.,  it  can  be  used 
;  evaluate  the  effect  of  design  parameter  changes  on  canceler  performance. 


14  SUBJECT  TERMS 

Adaptive  filter;  adaptive  noise  canceler;  cosite  interference 
reduction;  LMS  algorithm;  frequency  hopping  radio 


15.  NUMBER  OF  PAGES 

84 _ 

It.  PRICE  CODE 


,  *7  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

1  OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

j  Unclassified 

Unclassified 

Unclassified 

UL 

NSN  2S40  01  280  5S0O 


Standard  Form  298  (Rev  2-89) 

Pfpv:nb»d  by  ANSi  Ma  l 
2M-I02 


TABLE  OF  CONTENTS 


Section  Page 

INTRODUCTION  .  1 

ADAPTIVE  NOISE  CANCELING  .  3 

ADAPTIVE  NOISE  CANCELING  WITH  A  SINGLE  INPUT  .  8 

LMS  ALGORITHM .  11 

DESCRIPTION  OF  THE  ADAPTIVE  NOISE  CANCELER  WITH 
SINGLE  INPUT  SIMULATION  .  14 

SIMULATION  ASSUMPTIONS  .  16 

SIMULATION  INPUT  PARAMETERS  .  20 

DESIGN  EXAMPLE  .  28 

CONCLUSION .  75 

REFERENCES .  76 


iii 


A<JC933lon  1 

NT IS  GRA&I 
DTK:  TAB 

Un.-nuo-jj-iced 
Ju::t  if  icat ; 

por 

1ST 

a 

□ 

on 

Bv .  __ 

^Distribution/ 

_ Ayailabill 

>'  ■;  >  1  i 

foist  !  ,7  pr;  C 

if'1!  1 

ty  Codeo 

and/or 

ial 

LIST  OF  FIGURES 


Figure  Page 

1  Adaptive  Noise-Canceling  Concept  .  4 

2  Programmable  Transversal  Filter  (PTF)  .  5 

3  Adaptive  Noise  Canceler  with  a  Single  Input  .  9 

4  The  Adaptive  Linear  Combiner  .  13 

5  Simulation  Snapshot  Point  .  19 

6a  Iterated  Canceler  Output  Power  for  N  =  16  Taps, 

MU  =  0.1 .  35 

6b  PTF  Amplitude  Gain  vs  Frequency  for  N  -  16  Taps, 

MU  =  0  .  1 .  3  6 

6c  PTF  Phase  Shift  vs  Frequency  for  N  =  16  Taps, 

MU  =  0 . 1 .  37 

7a  Iterated  Canceler  Output  Power  for  N  =  32  Taps, 

MU  =  0.1 .  38 

7b  PTF  Amplitutde  Gain  vs  Frequency  for  N  =  32  Taps, 

MU  =  0. 1 .  39 

7c  PTF  Phase  Shift  vs  Frequency  for  N  =  32  Taps, 

MU  =  0.1 .  40 

3a  Iterated  Canceler  Output  Power  for  N  =  64  Taps, 

MU  =  0.1 .  41 

8b  PTF  Amplitude  Gain  vs  Frequency  for  N  =  64  Taps, 

MU  =  0.1 .  42 

3c  PTF  Phase  Shift  vs  Frequency  for  N  =  64  Taps, 

MU  =  0 . 1 .  43 

9a  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0. 1 .  44 

9b  PTF  Amplitude  Gain  vs  Frequency  for  N  =  128  Taps, 

MU  =  0. 1 .  45 

9c  PTF  Phase  Shift  vs  Frequency  for  N  =  1?.<J  Taps, 

MU  =  0.1 . 46 

10a  Iterated  Canceler  Output  Power  for  N  =  256  Taps, 

MU  =  0.04 .  47 


IV 


10b  PTF  Amplitude  Gain  vs  Frequency  for  N  =  256  Taps, 


MU  =  0.04 .  48 

10c  PTF  Phase  shift  vs  Freauency  for  N  =  256  Taps, 

MU  =  0.04 .  49 

11  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.2  55 

12  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.16 .  57 

13  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.15 .  59 

14  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.09 .  61 

15  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.08 .  62 

16  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.775  .  63 


17  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  -  0.075  .  64 

18  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.07 .  65 

19  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.06 .  66 

20a  PTF  Amplitude  Gain  vs  Frequency  for  N  =  128  Taps, 

MU  =  0.03 .  67 

20b  PTF  Phase  Shift  vs  Frequency  for  N  =  128  Taps, 

MU  =  0.08 .  68 

21  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.08,  and  Interfering  Signal 

Strength  =  10  dBm .  70 

22  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 

MU  =  0.8,  and  Interfering  Signal 

Strength  =  10  dBm .  73 


v 


LIST  OF  TABLES 

Table  Page 

1  Convergence  Parameter  vs  Number  of  Iterations 

Necessary  for  a  30  db  Reduction  in  Interference 

Power .  53 


vi 


INTRODUCTION 


This  report  presents  the  results  of  a  computer  simulation  of 
the  operation  of  an  adaptive  noise  canceler  in  the  simultaneous 
presence  of  both  strong  and  weak  narrow-band  signals.  The 
simulation  is  intended  to  be  used  as  a  tool  for  designing  such  a 
canceler.  For  a  given  scenario,  it  will  identify  the  design 
parameter  values  needed  to  cause  the  canceler  to  attenuate  the 
strong  signal  to  a  user  specified  level  and  pass  the  weak  signal. 

This  study  was  motivated  by  the  cosite  interference  problem 
encountered  by  frequency  hopping  radios.  When  two  or  more  such 
radios  and  their  antennas  are  independently  operated  in  close 
proximity,  i.e.,  in  a  jeep  or  communication  shelter,  a  cosite 
interference  problem  can  develop.  In  this  type  of  situation,  the 
radio  may  not  be  able  to  meet  its  specified  bit  error  rate.  A 
degraded  bit  error  rate  means  that  the  receiver  sensitivity  will 
be  degraded  and,  as  a  result,  communication  range  will  be  de¬ 
creased  . 

This  type  of  interference  problem  is  caused  by  the  trans¬ 
mitter's  strong  signal  being  too  close  to  the  frequency  the 
receiver  is  tuned  to.  The  difference  in  power  levels  between  the 
interfering  transmitter  signal  at  the  receiver  input  and  the 
"".nimuiri  signal  the  receiver  is  capable  of  detecting  could  be  in 
excess  of  130  dB.  The  receiver  may  not  be  able  to  provide  the 
entire  130  dB  of  interference  rejection  filtering  needed  at  the 
transmitter  frequency.  Therefore,  an  external  applique  capable 
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of  supplying  the  additional  filtering  may  be  required.  An  adap¬ 
tive  noise  canceler  with  a  single  input  is  one  possible  way  of 
providing  this  required  additional  filtering.  The  next  two 
sections  describe  the  adaptive  noise  canceling  concept. 
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ADAPTIVE  NOISE  CANCELING 


An  Adaptive  Noise  Canceler  is  shown  in  Figure  1.  It  works 
as  follows. 

"A  signal  is  transmitted  over  a  channel  to  a  sensor 
that  receives  the  signal  plus  an  uncorrelated  noise  NQ. 

The  combined  signal  and  noise  S+NQ  form  the  'primary 
input'  to  the  canceler.  A  second  sensor  receives  a 
noise  N ^  which  is  uncorrelated  with  the  signal  but 
correlated  in  some  unknown  way  with  the  noise  NQ.  This 
sensor  provides  the  'reference  input'  to  the  canceler. 

The  noise  is  filtered  to  produce  an  output  Y  that  is 
a  close  replica  of  NQ.  This  output  is  subtracted  from 
the  primary  input  S+N0  to  produce  the  system  output 
Nq  -  Y".1 

The  system  output  E  =  S  +  NQ  -  Y  of  the  canceler  is  used  to 
modify,  via  an  adaptive  algorithm,  the  frequency  response  of  the 
adaptive  filter. 

The  adaptive  filter  will  usually  be  implemented  as  a  pro¬ 
grammable  transversal  filter  (PTF)  (see  Figure  2).  A  transversal 
filter  is  the  preferred  implementation  because: 

1.  It  is  one  of  the  simplest  filter  structures.  The 
filter  output  is  simply  the  sum  of  delayed  and  scaled 
inputs . 

2.  There  is  no  feedback  from  the  taps  to  the  input. 

3.  It  is  stable.  Since  there  is  no  feedback,  a  finite 
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Figure  1.  Adaptive  Noise-Canceling  Concept  (from  Ref.  9). 


filter  input  produces  a  finite  filter  output.  If 
feedback  were  present,  then  an  inappropriate  choice  of 
filter  tap  weights  could  cause  the  filter  output  to 
become  unbounded,  i.e.,  the  output  would  oscillate.  A 
transversal  filter  cannot  oscillate. 

4.  It  has  a  linear  phase  characterist ic ,  i.e.,  it  produces  a 
phase  shift  that  is  linearly  proportional  to  frequency.  It 
can  be  shown14  that  if  a  signal  is  to  be  passed  through  a 
linear  system  without  any  resultant  distortion,  the  overall 
system  frequency  response  must  have  a  constant  amplitude 
gain  characteristic  over  the  frequency  spectrum  of  the  input 
signal  and  its  phase  shift  must  be  linear  over  the  same 
frequency  spectrum.  Filtering  without  distortion  is  impor¬ 
tant  for  adaptive  noise  canceling  because  the  adaptive 
filter  must  pass  the  interferer  without  distortion  so  that 
it  can  be  subtracted  (at  the  summer)  from  the  unfiltered 
interferer.  If  the  adaptive  filter  introduces  distortion, 
then  the  summer  is  no  longer  subtracting  two  identical 
interferers.  As  a  result,  the  error  power  and  error  ampli¬ 
tude  with  adaptive  filter  distortion  will  be  higher  than 
without  distortion.  This  is  not  desirable. 

5.  There  is  a  simple  and  analytically  tractable  relation¬ 
ship  between  the  transfer  function  of  a  transversal 
filter  and  its  parameters  (see  Equation  2  of  the  simu¬ 
lation  description  section).  The  complicated  nonlinear 
relationship  between  parameters  and  transfer  function 
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for  most  other  filter  structures  makes  the  analysis  and 
calculation  of  adaptive  algorithms  much  more  difficult 
than  for  transversal  filters. 

6.  Widrow's  algorithm,  one  of  the  most  widely  used  adaptive 
algorithms,  assumes  a  transversal  filter  structure. 

A  PTF  forms  a  weighted  sum  of  delayed  versions  of  the  input 
signal.  It  is  programmable  in  that  the  weights  can  be  changed. 
Changing  the  weights  changes  the  frequency  transfer  function  of 
the  PTF.  A  PTF  is  identical  in  structure  to  a  programmable 
finite  impulses  response  (FIR)  digital  filter.  The  specific 
technology  used  to  implement  a  PTF  will  depend  on  the  frequency 
range  of  interest.  For  VHF  and  UHF  applications,  Surface  Acous¬ 
tic  Wave  (SAW)  devices  are  an  appropriate  technology.  At  VHF  and 
UHF  frequencies,  SAW  technology  can  give  the  appropriate  sam¬ 
pling  rates  (intertap  delay)  and  total  delay  times  necessary  to 
implement  transversal  filters  with  the  required  bandwidth  and 
frequency  resolution  needed  for  cosite  interference  reduction. 
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ADAPTIVE  NOISE  CANCELING  WITH  A  SINGLE  INPUT 


Before  an  adaptive  noise  canceler  can  be  implemented,  a 
reference  signal  correlated  with  the  interfering  signal  but  not 
the  intended  signal  must  be  generated.  When  the  interfering 
signal  NQ  is  much  stronger  than  the  intended  signal  S,  the 
reference  signal  can  be  generated  by  modifying  the  adaptive  noise 
canceler  of  Figure  1  to  give  the  circuit  shown  in  Figure  3  (an 
adaptive  noise  canceler  with  a  single  input).  In  Figure  3,  the 
primary  and  reference  inputs  are  connected  together.  In  effect, 
Figure  3  assumes  that  the  reference  input  is  equal  to  the  primary 
input.  This  may  at  first  appear  contradictory.  The  reference 
input  Nj_  (see  Figure  1)  has  to  be  correlated  to  the  interference 
N0,  not  the  signal  S.  But  since  the  signal  S  is  part  of  the 
primary  input,  it  will  be  part  of  reference  input  if  reference 
input  =  primary  input  as  per  Figure  3.  Hence  the  reference  input 
appears  to  be  correlated  to  the  signal,  also. 

When  the  interfering  signal  N0  is  much  larger  than  the 
intended  signal  (NQ>>S)  the  apparent  contradiction  is  resolved. 

In  this  case  the  reference  input  (N^  =  S  +  N0  =  primary 

input)  is  highly  correlated  with  and  "looks"  like  the  interfering 
signal  (i.e.,  ~  N0)  .  While  S  is  a  component  of  N ^ ,  and 

therefore,  will  correlate  to  a  certain  extent  with  N-j_,  NQ  is  so 
much  larger  than  S  that  will  be  much  more  highly  correlated  to 
Nq  than  S.  So  to  a  very  good  approximation,  the  reference  input 
N]  is  correlated  to  the  interference  NQ,  not  the  signal  S.  (For 
a  rigorous  proof  ,  Reference  2.) 
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Figure  3.  Adaptive  Noise  Canceler  with  a  Single  Input. 


It  will  now  be  shown  why  the  reference  input  must  be  corre¬ 
lated  to  the  interference  and  not  the  signal.  The  adaptive 
filter  within  the  canceler  must  filter  the  reference  input  to 
produce  an  output  Y  that  is  a  close  replica  of  N0.  If  is  not 
correlated  to  NQ,  i.e.,  if  N]_  does  not  "look"  somewhat  like  NQ, 
then  no  amount  of  filtering  can  make  Y  look  like  N0. 

The  factor  1//2  appears  in  Figure  3  because  the  input  power 
splitter  is  assumed  to  evenly  split  the  power  associated  with  the 
signal  and  interference  amplitudes  S  and  NQ .  Since  power  is 
proportional  to  amplitude  squared,  reducing  power  by  a  factor  of 
two  means  that  amplitude  is  reduced  by  J2  at  each  output  of  the 
input  power  splitter. 

As  the  adaptive  algorithm  iterates,  it  will  cause  the 
adaptive  filter  to  form  a  bandpass  around  the  interfering  fre¬ 
quency.  If  the  PTF  has  been  properly  designed,  then  the  result¬ 
ing  bandpass  filter  will  "pass"  FNq  the  interfering  frequency  and 
"reject"  the  intended  signal  frequency.  Then  the  output  of  the 
adaptive  filter  (the  filtered  reference  signal)  will  "look"  even 
mere  like  NQ  than  input  signal.  When  the  output  is  subtracted 
from  S  +  Nq,  at  the  output  power  combiner,  a  signal  very  similar 
to  S  will  remain.  The  interference  has  been  canceled.  Thus  the 
circuit  shown  in  Figure  3  does  indeed  behave  as  an  adaptive  noise 
canceler. 
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LMS  ALGORITHM 


Adaptive  algorithms  have  been  extensively  discussed  in  a 
previous  LABCOM  technical  report2,  in  books-*-  >  3  >  ^  >  6  <  7  >  ®  and 
review  articles.9'-1-0  The  theory  of  adaptive  algorithms  is  too 
lengthy  and  involved  to  be  reviewed  in  this  report.  The  inter¬ 
ested  reader  is  referred  to  the  references  just  cited. 

The  Least  Mean  Square  (LMS)  or  Widrow's  algorithm  was  chosen 
as  the  adaptive  algorithm  for  use  in  the  simulation.  It  is  well 
understood,  computationally  simple  and  fast. 

In  the  LMS  algorithm,  it  is  assumed  that  the  adaptive  filter 
is  an  adaptive  linear  combiner  (see  Figure  4).  If  data  are 
acquired  and  input  in  parallel  to  an  adaptive  linear  combiner, 
the  structure  in  Figure  4a  is  used.  For  serial  data  input,  the 
structure  in  4b  is  used.  Note  that  Figure  4b  is  just  a  tapped 
delay  line  or  transversal  filter.  This  is  the  filter  structure 
that  will  be  assumed  in  this  report.  It  is  further  assumed  that 
a  "desired"  response  signal  is  available.  For  the  adaptive  noise 
cancoler  with  a  single  input  of  Figure  3,  the  primary  input  is 
the  "desired"  response. 

The  LMS  algorithm  is  given  by: 

wk+l  =  wk  +  2^ekxk  (!) 

where : 

wK  =  the  weight  vector  at  the  kth  iteration,  i.e.,  the  set  of 
tap  weights  used  on  the  kth  iteration. 
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Wkfl  =  the  weight  vector  at  the  k+lth  iteration. 

M  =  a  constant  that  regulates  the  step  size  or  increment  size 
of  the  weight  vector  change. 

=  canceler  output  error  amplitude  at  the  kth  iteration. 

Xk  =  tap  signal  vector  at  the  kth  iteration,  i.e.,  the  set  of  tap 
signal  amplitudes  at  the  kth  iteration. 
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Figure  4.  The  Adaptive  Linear  Combiner:  (a)  In  General  Form; 
(b)  As  a  Transversal  Filter  (from  Ref.  1 ,  page  101). 


DESCRIPTION  OF  THE  ADAPTIVE  NOISE  CANCELER  WITH 
SINGLE  INPUT  SIMULATION 


The  circuit  in  Figure  3  was  simulated.  In  the  simulation, 

two  sinusoids  of  different  frequency  (user  specified)  were  input 

to  the  power  splitter.  For  any  given  iteration  k,  the  weight 

vector  for  that  iteration  was  used  to  calculate  the  frequency 

transfer  function  of  the  adaptive  filter  or  PTF.  Knowledge  of 

the  PTF  frequency  transfer  function  and  the  input  frequencies 

allows  the  input  to  the  negative  or  inverting  part  of  the  power 

combiner  from  the  lower  "channel"  of  the  adaptive  noise  canceler 

to  be  calculated.  The  output  of  the  power  combiner,  the  error 

amplitude  and  the  error  power  Pe  can  then  be  calculated.  The 

signals  Xj,  at  the  tap  outputs  of  the  PTF  are  then  calculated. 

Since  and  are  known,  Widrow's  LMS  algorithms  (Equation  1) 

can  be  used  to  calculate  a  new  set  of  tap  weights  Wjc+1.  The 

simulation  can  then  be  repeated  with  the  new  tap  weight  vector 
— * 

Wk+i  and  new  values  for  the  error  amplitude  and  the  error 

power  PPl  can  then  be  calculated.  After  each  iteration,  the 
c  k 

adaptive  noise  canceler  error  or  output  power  for  that  iteration 
is  displayed.  When  the  output  power  has  been  reduced  by  a  user 
specified  number  of  dB,  the  simulation  is  stopped. 

Once  the  adaptive  filter  processor  speed  of  operation  or 
equivalently  its  tine  per  iteration  is  specified,  then  the 
product  of  the  number  of  iterations  necessary  to  reduce  the 
output  power  by  the  user  specified  number  of  dB  and  the  time  per 
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iteration  gives  the  adaptive  convergence  time.  This  is  the  time 
needed  for  the  required  reduction  in  output  power  to  take  place 
So,  in  effect,  when  the  simulation  stops,  the  number  of  itera¬ 
tions  it  has  performed  tells  us  the  value  of  the  adaptive  conver 
gence  time. 

The  simulation  also  displays  the  PTF  frequency  transfer 
function.  In  order  for  the  canceler  to  work  properly  the  inter¬ 
fering  frequency  must  be  "passed"  and  the  intended  signal  "re¬ 
jected"  by  the  PTF.  This  can  be  observed  directly  from  the  PTF 
amplitude  gain  vs.  frequency  curve. 
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SIMULATION  ASSUMPTIONS 


The  assumptions  made  in  the  course  of  the  simulation  are  as 
follows : 

1.  The  programmable  transversal  filter  was  assumed  to  have 
equally  spaced  taps.  The  frequency  transfer  function  for  a 
transversal  filter  with  uniformly  spaced  taps  is  given  by 
the  following  equations: 

N 

Hk(o>)  = .  Z  Wik  e3“A(-i)  (2) 

where : 

Hk(o)  =  frequency  transfer  function  at  the  kth  iteration. 

w  =  frequency 

wik  =  weight  on  the  ith  tap  at  the  kth  iteration. 

A  =  intertap  delay 

N  =  number  of  taps 

j  =  f= 1 

amplitude  gain  =  jHk(o)j  =  [(Real  Hk(o)))^  +  (Imag  Hk  (u)  ) 2  ]  V  2  (3) 
phase  shift  =  tan-1  (Imag  Hk((j)  /  Real  Hk(w))  (4) 

2.  The  error  amplitude  output  from  the  adaptive  noise 
canceler  was  calculated  using  a  "snapshot  in  time" 
approach,  i.e.,  the  interfering  and  intended  signals 
input  to  the  power  splitter  of  the  adaptive  noise 
canceler  were  calculated  at  a  single  point  in  time. 
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These  signals  are  then  mathematically  sent  through  the 
adaptive  noise  canceler.  The  upper  channel  of  the  canceler  does 
not  change  either  sinusoid.  The  PTF  in  the  lower  channel  changes 
the  amplitude  and  phase  of  each  signal.  The  two  channels  are 
then  subtracted  at  the  power  combiner  to  produce  the  error 
amplitude 

3.  At  the  input  to  the  power  splitter  the  intended  and  inter¬ 
fering  signals  are  assumed  to  be  in  phase.  An  input  parame¬ 
ter  (called  SIMTIME)  determines  at  what  point  in  time  the 
snapshot  is  taken  (see  Figure  5). 

4.  At  each  iteration  it  was  assumed  that  the  snapshot  was  taken 
at  the  same  relative  point  on  the  waveform,  i.e.,  the  same 
value  of  SIMTIME  was  used  at  each  iteration. 

5.  It  was  found  that  the  only  way  that  the  canceler  output 
could  be  reduced  by  the  specified  number  of  dB  was  to  choose 
the  snapshot  point  very  close  (within  a  few  degrees)  to  the 
maximum  value  of  the  interference.  The  snapshot  point  for 
all  the  simulation  runs  was  therefore  chosen  at  the  maximum 
value  of  the  interference.  This  is  not  as  restrictive  an 
assumption  as  it  appears  to  be  and  will  be  further  discussed 
in  the  parameters  section  of  this  report. 

o.  The  output  power  for  the  power  combiner  was  calculated  as 
the  RMS  output  of  the  two  inputs. 

7.  The  first  or  initial  set  of  tap  weights  is  always  assumed  to 
be  all  zero's.  This  was  done  so  that  all  simulation  runs 
have  the  same  "starting  point,"  i.e.,  initial  weight  vector 
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(0,  0...,  0).  Thus,  meaningful  comparisons  between  runs  can 
be  made . 

8.  The  simulation  treats  the  two  frequency  hopping  radios  as 

fixed  frequency  radios,  i.e.,  it  does  not  hop  or  change  the 
frequencies  input  by  the  user.  This  assumption  was  made  in 
order  to  observe  how  many  iterations  of  the  adaptive  algo¬ 
rithm  were  needed  to  reduce  the  interfering  signal  by  the 
user  specified  number  of  dB.  If  the  radios'  frequencies 
were  allowed  to  hop,  then  for  a  "bad"  choice  of  design 
parameters  the  hop  would  occur  before  the  power  had  been 
reduced  by  the  specified  amount.  It  would  not  be  possible 
to  tell  just  how  bad  the  design  parameters  were.  Frequency 
hopping  would  only  complicate  the  simulation,  it  would  not 
help  in  understanding  adaptive  algorithm  perf ormance . 
Therefore,  it  was  not  simulated. 
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Figure  5.  Simulation  Snapshot  Point. 


SIMULATION  INPUT  PARAMETERS 


The  user  supplied  input  parameters  to  the  simulation  can  be 
divided  into  two  categories:  scenario  parameters  and  design 
parameters . 

Scenario  parameters  are  fixed  or  given  by  the  exact  nature 
of  the  scenario  simulated.  The  analyst  has  no  control  over 
these  parameters.  Specifically  the  scenario  input  parameters 
required  are: 

o  Intended  frequency  and  power. 

o  Interfering  frequency  and  power. 

o  Desired  reduction  of  interfering  signal  power  (in  dB) . 

o  PTF  voltage  loss  per  tap. 

Design  parameters  can  be  varied.  They  are  chosen  to  achieve 
the  required  adaptive  noise  canceler  frequency  transfer  function 
and  adaptive  convergence  time. 

The  design  input  parameters  required  are: 

o  The  intertap  delay  time  between  taps  of  the  PTF. 

o  The  number  of  taps  in  the  PTF. 

o  MU  (or  11)  the  convergence  parameter. 

A  The  intertap  delay  A  should  be  chosen  to  be  consistent 
with  the  Nyquist  criteria,  i.e.,  if  fmax  is  the  highest  frequency 
in  the  baseband  signal,  then 

f max  -  — iL-  (5) 

2  A 
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N,  the  number  of  taps,  should  be  chosen  on  the  basis  of  the 
required  PTF  frequency  resolution,  where  frequency  resolution  is 
defined  as  the  difference  in  frequency  between  the  center  of  the 
PTF  frequency  transfer  function  and  its  first  zero.  Frequency 
resolution  is  given  by11: 

Frequency  Resolution  =  1/NA  =  1/ (total  delay)  (6) 

The  convergence  parameter  jit,  which  regulates  the  size  of  the 
weight  vector  increment  at  each  iteration  (see  Equation  1)  is 
the  most  important  design  parameter.  Adaptive  convergence  time 
is  a  very  sensitive  function  of  p.  If  M  is  too  small,  the 
convergence  time  will  be  too  long.  If  n  is  too  large,  the  LMS 
algorithm  may  "blow  up"  and  no  convergence  will  be  achieved, 
i.e.,  there  will  be  no  reduction  in  the  interfering  signal  power. 

It  can  be  shown12  that  for  a  transversal  filter  with  no  tap 
loss,  if  the  following  inequality  is  satisfied: 

o<  m  <  1  /  ((N+l)  (input  power))  (7) 

then  the  expected  adaptive  noise  canceler  weight  vector  will 
converge  to  the  optimum  or  Weiner  weight  vector.  The  simulation 
requirements  are  somewhat  different.  Total  elimination  of  the 
interference  (the  Wiener  solution)  is  not  required.  It  is  only 
necessary  to  reduce  the  interfering  signal  by  a  user  specified 
iiumoer  of  dB,  but  it  must  be  done  in  a  minimal  number  of  itera¬ 
tions.  Despite  these  differences,  Inequality  7  has  been  found 
(empirically)  to  be  a  good  guide  for  the  selection  of  n . 

If  the  particular  technology  used  has  a  loss  at  the  PTF  taps 
(e.g.,  SAW  technology)  then  n  as  determined  from  Inequality  7 
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should  be  multiplied  by  (1/tap  loss)  to  compensate  for  the  tap 
loss.  This  works  because  in  the  LMS  algorithm  (Equation  1) 


wk+l  ~  wk  +  2^ekxk  f1) 

— ♦ 

H  and  the  tap  signal  vector  only  appear  in  a  product.  If 
— * 

Xk  is  reduced  by  the  same  loss  at  each  tap  then  multiplying  n  by 
(1/tap  loss)  cancels  the  tap  loss  in  the  product.  So  with  no 
loss  in  generality  in  the  analysis,  it  can  be  assumed  that  tap 
loss  =  0.  Therefore  the  analysis  can  be  performed  before  tap 
loss  data  is  available. 

The  only  remaining  parameter  is  SIMTIME,  the  parameter  that 
determines  the  point  in  time  at  which  the  simulation  snapshot  is 
taken  (see  the  simulation  assumption  section  and  Figure  5) . 

Using  typical  cosite  interference  input  parameters  it  was  found 
that  the  only  way  that  the  canceler  output  could  be  reduced  30  dB 
(a  "typical"  cosite  interference  requirement)  was  to  choose  the 
snapshot  point  very  close  (within  a  few  degrees)  to  the  maximum 
amplitude  value  of  the  interference.  So,  in  effect,  SIMTIME  is 
no  longer  an  adjustable  input  parameter.  All  simulation  runs 
will  use  the  same  value  (specifically,  SIMTIME  =  1) . 

Choosing  the  snapshot  point  at  the  maximum  amplitude  value 
of  Ihe  interfering  signal  can  be  justified  on  both  mathematical 
and  physical  grounds. 

Mathematically,  the  dominant  contribution  to  canceler  output 
is  from  the  interfering  signal.  Any  signal  that  is  input  to  an 
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adaptive  noise  canceler  with  a  single  input  (see  Figure  3)  is 
first  split  by  the  power  splitter.  Half  the  input  power  goes 
directly  to  the  non-inverting  input  of  the  power  combiner.  The 
other  half  of  the  power  is  first  filtered  by  the  PTF  and  then 
input  to  the  inverting  input  of  the  power  combiner.  In  the 
course  of  the  simulation,  the  canceler  output  power  due  to  the 
interfering  signal  was  shown  to  be  equal  to  and  was  modeled  as: 
Canceler  output  power  due  to  interfering  signal 

=  [c2J  +  C2j  |Hk(G>j)|2  -  2  C2j  |Hk(oj)|  Cos  (PSJSFk)]/4  (8) 

=  C2j  [1  +  | Hk (oj) |  (|Hk(oj)|  -  2Cos  (PSJSFk))]  /  4 

where : 

Cj  =  [2  •  Interfering  power  in  watts] ^ 

Hk(oj)  =  frequency  transfer  function  at  the  kth  iteration 

evaluated  at  the  interfering  or  jamming  frequency  uj. 

PSJSFk  =  PTF  phase  shift  at  the  jammer  signal  frequency  at  the 
kth  iteration. 

Not  being  able  to  achieve  the  user  specified  reduction  in 
the  interfering  power  (when  the  simulation  snapshot  is  not  taken 
at  the  interference  maximum)  means  that  the  power  calculated  in 
Equation  8  is  too  large.  This  means  that  either  the  PTF  phase 
shift  at  the  interfering  frequency  is  too  large  or  that  |Hk(oj)| 
i  5  not  numerically  close  enough  to  one.  |Hk(o>j)|,  the  PTF 
orplitude  gain,  must  be  close  to  one  and  the  phase  shift  must  be 
close  to  zero  in  order  that  the  filtered  interfering  frequency 
signal  cancels  the  unfiltered  interfering  frequency  at  the  power 
combiner . 
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The  PTF  phase  shift  is  given  by  Equation  4 


PSJSFk  =  tan-1  (Imag  Hk(wj)  /  Real  Hk(wj) )  (4) 

where  Hk(oj),  the  frequency  transfer  function  at  the  interfering 
frequency  oj ,  is  given  by  Equation  2. 

Hk  («j)  =JX  wI,k  (2) 

where : 

WI  k  “  Weight  on  the  Ihh  tap  at  the  kth  iteration 
(jj  =  Interfering  frequency 
A  =  Intertap  delay 
N  =  Number  of  taps 

j  =  y^r 

In  Equations  4  and  2,  both  the  PTF  phase  shift  and  amplitude 
gain  are  functions  of  the  tap  weights.  The  tap  weights  are  given 
by  Widrow's  LdS  algorithm. 

— *  — ► 

wk+l  =  wk  +  2^ekxk  (!) 

The  value  of  the  error  amplitude  ek  will  depend  on  the 
snapshot  time,  i.e.,  ck  will  be  a  function  of  SIMTIME.  ek  is  the 
amplitude  output  of  the  power  combiner.  Since  the  interfering 
signal  is  so  much  larger  than  the  intended  signal,  ek  can  be 
approximated  by  the  difference  between  the  unfiltered  and  fil¬ 
tered  interfering  signals.  There  is  a  phase  shift  between  the 
unfiltered  and  filtered  interfering  signals  that  is  given  by 
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Equation  4.  This  phase  shift  affects  the  value  of  ek- 

Choosing  the  snapshot  point  at  the  maximum  value  of  the 
interfering  signal  (i.e.,  SIMTIME  =  1)  minimizes  the  impact  of  a 
nonzero  phase  shift  on  ek-  In  other  words,  when  SIMTIME  =  1,  we 
can  achieve  a  "small"  value  of  ek  even  for  a  non-zero  phase 
shift.  This  is  due  to  the  fact  that  at  the  maximum  value  of  the 
interference,  i.e.,  at  the  top  of  the  sinewave,  the  sinewave  is 
flat.  Its  first  derivative  is  zero.  Therefore,  its  rate  of 
change  is  also  zero.  This  is  the  point  or  region  where  the 
sinewave  changes  most  slowly.  This  region  of  slow  change  is 
where  a  "large"  phase  shift  can  be  most  easily  "tolerated,"  i.e., 
a  phase  shift  produces  the  least  change  in  the  sinewave  when  the 
snapshot  point  is  at  the  interfering  frequency  maximum.  The 
phase  shift  does  not  move  the  sinewave  very  far  off  of  its 
maximum.  This  implies  that  €k  will  be  small.  So,  in  effect, 
SIMTIME  =  1  produces  the  smallest  error  amplitude  ek  for  any 
given  phase  shift. 

If  ek  is  "small"  (when  SIMTIME  =  1),  then  in  Widrow's 
a  Igor i thm 

Wk+i  -  Wk  +  2^ekXk  (1) 

— ♦ 

cue  incremented  term  2/i  ek  Xk  will  also  be  relatively  small.  A 
smaller  incremental  term  means  that  the  optimal  set  of  weights 
can  be  approached  more  precisely,  i.e.,  Equation  1  can  get  closer 
co  the  optimal  set  of  weights  with  a  "small"  incremental  term 
than  with  a  "large"  term.  As  we  approach  the  optimal  weight 
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vector,  |Hk(«j) |  approaches  closer  to  one  (via  Equation  2)  and 
PSJSFk  approaches  closer  to  zero  (via  Equation  4).  This  will  in 
turn  reduce  (via  Equation  8)  the  adaptive  noise  canceler  output 
power  due  to  the  interfering  signal.  This  is  why  the  simulation 
worked  best  when  SIMTIME  =  1.  To  summarize:  choosing  the 
snapshot  point  at  the  interfering  signal  maximum  (SIMTIME  =  1) 
gives  a  "small"  error  output  ck.  This  gives  a  "small"  incremen¬ 
tal  term  2 uckXk  in  Widrow's  algorithm.  This  small  incremental 
term  lets  Equation  1  get  "very  close"  to  the  optimal  tap  weight 
vector.  As  we  get  close  to  the  optimal  tap  weight  vector, 

|Hk(uj)  I  1  and  PSJFSk  -  0.  These  two  conditions  reduce  the 
canceler  output  power  to  the  extent  necessary  to  achieve  the  user 
specified  power  reduction. 

Choosing  the  snapshot  point  at  the  maximum  value  of  the 
interfering  signal  can  also  be  justified  on  physical  grounds. 

Real  signal  measurements  are  not  made  at  a  single  point  in  time 
(as  assumed  in  the  simulation) .  They  are  made  over  an  interval 
of  time.  Amplitude  or  power  measurements  are  integrated  over 
this  period  of  time  by  the  measuring  instrument.  This  integra¬ 
tion  insures  that  measurements  of  strong  signals  are  much  larger 
than  measurements  of  weak  signals.  Any  simulation  of  a  real 
measurement  must  be  able  to  reproduce  this  behavior.  Choosing 
the  snapshot  time  in  the  simulation  at  the  maximum  value  of  the 
interfering  sinewave  guarantees  this  behavior.  This  is  so 
because  in  the  scenario  the  interfering  signal  will  be  much 
larger  than  the  intended  signal.  This  implies  that  the  maximum 
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value  of  the  strong  amplitude  will  be  much  larger  than  the 
amplitude  of  the  weak  intended  signal  at  the  snapshot  time.  This 
is  the  behavior  the  simulation  must  and  does  reproduce. 

If  the  snapshot  time  did  not  correspond  to  the  maximum  value 
of  the  interferer,  then  the  dynamic  range  between  the  interferer 
and  the  intended  signal  would  be  compressed.  If,  for  example,  a 
randomly  chosen  snapshot  point  (actually  a  snapshot  time)  on  the 
interferer  is  used  instead  of  using  the  maximum  value,  then  the 
"difference"  between  interfering  and  intended  signals  at  the 
random  point  would  not  be  as  large  as  the  "difference"  calculated 
by  assuming  that  the  snapshot  was  taken  at  the  interference 
maximum.  Choosing  the  snapshot  time  at  the  maximum  value  of  the 
interference  is  a  way  of  insuring  that  the  simulation  can 
distinguish  between  strong  and  weak  signals  without  compressing 
the  input  signal  dynamic  range. 

As  further  physical  justification,  consider  that  any  real 
signal  measurement,  whether  defined  via  an  average  or  RMS  proce¬ 
dure,  will  be  proportional  to  the  maximum  value  of  the  waveform; 
i.e.,  for  a  sinusoid  Vpsin(oT)  (where  Vp  is  the  maximum  value  of 
the  waveform)  the  average  value  is  (2/7r)Vp,  the  RMS  value  is 
Vp//2.  So  by  choosing  in  the  simulation  the  maximum  value  (Vp) 
of  me  interfering  signal  as  our  snapshot  point  we  are  getting  an 
insight  into  the  real  signal  behavior. 
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DESIGN  EXAMPLE 


In  order  to  illustrate  the  use  of  the  adaptive  noise  cancel- 
er  simulation  it  will  be  applied  to  the  analysis  of  two  cosited 
frequency  hoppinq  radios.  The  canceler  will  be  used  to  protect 
the  receiving  radio  from  the  transmitting  radio's  signal.  The 
simulation  will  be  used  to  arrive  at  a  set  of  design  parameters 
for  the  canceler  that  minimizes  the  number  of  iterations  neces¬ 
sary  to  achieve  the  user  specified  power  reduction  in  the  trans¬ 
mitter  or  interfering  frequency  signal  while  at  the  same  time 
passing  the  receiver  intended  frequency. 

The  simulation  treats  the  two  frequency  hopping  radios  as 
fixed  radios.  It  will  be  assumed  that  the  transmitter  frequency 
differs  from  the  receiver  frequency  by  1  MHz  (i.e.,  the  user 
selected  intended  and  interfering  frequencies  will  be  chosen  1 
MHz  apart  and  then  input  to  the  simulation,  the  simulation  does 
not  assume  a  1  MHz  separation.)  A  one  MHz  frequency  difference 
was  chosen  because  a  significant  reduction  in  the  interfering 
signal  power  1  to  5  MHz  away  from  the  receiver  frequency  will 
significantly  reduce  the  receiver  bit  error  rate  (BER)  and 
increase  the  communications  range.  Transmitter  signals  less  than 
3  MHz  away  from  the  receiver  frequency  will  not  be  considered 
because  if  the  frequency  hopping  transmitter  is  hopping  randomly, 
these  "close-in"  frequencies  will  occur  so  infrequently  (compared 
to  transmitter  signal  frequencies  1-5  MHz  away  from  the  receiver 
frequency)  that  their  effect  on  the  receiver  BER  will  be  minimal. 
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In  a  typical  cosite  scenario  the  transmitter  power  will  be 
50  watts  (+47  dBm) .  The  required  receiver  sensitivity  will  be 
-98  dBm.  Propagation  loss  between  the  two  radio  antennas  will  be 
frequency  dependent.  It  can  vary  between  13  and  27  dB.  The 
canceler  is  intended  to  help  the  frequency  hopping  receiver 
filter  out  the  interference.  It  is  not  expected  to  do  all  the 
filtering  by  itself,  i.e.,  it  is  not  expected  to  reduce  the 
interferer  below  the  radio  sensitivity.  This  would  require  the 
canceler  to  produce  118-134  dB  or  rejection.  The  combination  of 
the  canceler  plus  the  tuned  circuitry  of  the  hopping  radio  is 
expected  to  reduce  the  interferer  below  the  radio  sensitivity.  A 
realistic  requirement  for  adaptive  noise  canceler  filtering  is  30 
dB.  SAW  technology  can  easily  achieve  30  dB  of  filtering. 
Together,  canceler  filtering  of  30  dB,  the  propagation  loss 
between  antennas,  and  the  radio  tuned  circuitry  will  reduce  the 
transmitter  signal  below  the  receiver  sensitivity. 

The  frequency  hopping  radios  will  be  assumed  to  hop  between 
30  and  90  MHz.  For  reasonable  insertion  loss,  linear  phase  SAW 
devices  are  limited  to  33%  fractional  bandwidth.  The  center 
frequency  of  the  SAW  device  in  the  canceler  will  therefore  be: 

fc  =  Bandwidth  =  60  MHz  =  180  MHz  (9) 

Fractional  Bandwidth  0.33 

Since  180  MHz  is  the  center  frequency  of  the  SAW  device  in 
the  PTF,  It  was  chosen  as  the  intended  signal  frequency  in  the 
simulation  examples  to  be  presented  in  this  report.  It  is  of 
course  assumed  that  the  hopping  radio  signals  will  be  up- 


converted  to  this  frequency  range,  processed  and  then  down- 
converted  . 

The  interference  power  input  to  the  receiver  (transmitter 
power  minus  propagation  loss)  will  be  20-34  dBm  depending  on 
frequency.  This  is  too  much  power  for  the  PTF/SAW  device  in  the 
canceler.  The  maximum  power  that  can  be  safely  input  to  it  is 
+20  dBm.  One  possible  way  of  protecting  the  SAW  device  is  by  use 
of  a  frequency  selective  limiter  that  filters  out  high  power 
signals  and  passes  low  power  signals.  Such  a  limiter  can  be 
simulated  by  simply  setting  the  input  interfering  power  at  +20 
dBm  and  the  input  intended  signal  at  the  receiver  sensitivity  of 
-98  dBm. 

There  is  now  enough  information  to  specify  the  scenario 
input  parameters: 

o  The  intended  frequency  will  be  assumed  to  be  at  180 

MHz,  the  center  frequency  of  the  SAW  device.  It  could 
have  been  chosen  anywhere  between  150  and  210  MHz. 

o  The  intended  signal  power  will  be  assumed  to  be  at  -98 
dBm,  the  lowest  power  the  receiver  is  capable  of  de¬ 
tecting  . 

o  The  interfering  frequency  will  be  assumed  to  be  at  181 
MHz,  one  MHz  away  from  the  intended  signal  at  180  MHz. 
It  could  have  been  assumed  the  interferer  was  at  179 
MHz.  It  makes  no  difference  whether  181  or  179  MHz  is 
chosen . 
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The  interfering  signal  power  will  be  assumed  to  be  at 
+  20  dBm,  the  output  of  a  frequency  selective  limiter. 

o  The  desired  reduction  of  the  interfering  signal  power 
will  be  30  dB. 

o  The  PTF  voltage  loss  per  tap  will  be  assumed  to  be 

zero.  (See  "Simulation  Input  Parameters"  section  for 
rationale . ) 

With  the  scenario  parameters  now  given,  the  simulation  will 
be  used  to  optimize  two  of  the  three  design  parameters  (the 
number  of  taps  in  the  PTF  and  n  the  convergence  parameter.)  The 
third  design  parameter,  the  intertap  delay  time  between  taps, 
will  be  assumed  to  be  6.9444  nanoseconds.  This  corresponds  to  a 
sampling  frequency  of  144  MHz.  This  sampling  frequency  was 
chosen  in  order  to  model  a  SAW/PTF  currently  being  built  for  ETDL 
by  Texas  Instrument under  Contract  DAAL01-88-C-083 1 .  In  that 
effort  a  sampling  frequency  of  144  MHz  is  being  used.  144  MHz  is 
-lightly  larger  than  the  Nyquist  sampling  rate  (120  MHz)  neces¬ 
sary  to  sample  a  60  MHz  bandwidth  signal.  It  was  chosen  by  TI  to 
provide  some  protection  against  aliasing. 

To  be  more  specific:  discrete  sampling  of  an  analog  wave¬ 
form,  which  is  what  the  PTF  taps  do  (they  form  discrete  samples 
of  the  input  analog  waveform) ,  not  only  duplicates  the  input 
spectrum,  but  also  replicates  it  around  harmonics  of  the  sampling 
rate.  If  the  sampling  rate  (equal  to  l/(intertap  delay  time))  is 
luo  low,  the  replicated  spectrums  overlap  in  the  frequency 
domain.  Thus,  any  frequencies  higher  than  half  the  Nyquist  rate 


that  are  present  will  be  aliased,  i.e.,  appear  (due  to  undersamp- 
ling)  as  lower  frequencies.  This  distorts  or  equivalently 
introduces  interference  into  the  ••lower"  (below  half  the  Nyquist 
rate)  frequency  spectrum.  The  minimum  sampling  rate  necessary  to 
assure  no  overlap  in  the  output  spectrum  is  the  Nyquist  rate 
(equal  to  twice  the  highest  frequency  present  in  the  input) .  To 
insure  that  the  replicated  spectrums  do  not  overlap,  it  is 
considered  good  engineering  practice  to  sample  at  a  somewhat 
higher  rate  than  the  Nyquist  rate.  This  is  why  144  MHz  was 
chosen  as  the  sampling  frequency  rather  than  120  MHz. 

The  simulation  will  now  be  used  to  determine  N  the  number  of 
taps  needed  in  the  SAW/PTF.  Although  the  optimum  value  of  /x  the 
convergence  parameter  has  not  as  yet  been  determined,  a  value  to 
input  as  a  design  parameter  is  still  needed.  The  value  of  /x  that 
is  used  need  not  be  its  optimum  value  (that  will  be  determined 
after  the  number  of  taps  is  determined) .  The  convergence  parame¬ 
ter,  \i,  need  only  be  close  enough  to  the  optimum  (for  a  given  N) 
so  that  the  canceler  output  does  not  "blow  up",  but  eventually 
converges.  The  adaptive  convergence  time  is  affected  by  /x,  but  /i 
does  not  affect  the  "optimum"  frequency  transfer  function  of  the 
PTF .  For  the  purpose  of  determining  a  suitable  value  for  N  it 
does  not  matter  how  long  it  takes  to  arrive  at  the  optimum 
frequency  transfer  function.  Inequality  (7) 

1 

0  ...  n  <  -  (7) 

(N+l) (Input  Power) 
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can  be  used  as  a  guide  to  guessing  one  or  more  values  of  4  to  use 
in  determining  a  suitable  value  of  N.  For  the  given  scenario,  as 
N  varies  between  16  and  256,  4,  as  given  by  Inequality  7,  will 
vary  between  0.625  and  0.039.  Either  Inequality  7  can  be  used  to 
calculate  a  close  to  optimum  4  for  each  value  of  N  (equal  to  the 
upper  bound  of  the  inequality) ,  or  an  educated  guess  at  4  between 
0.625  and  0.039  can  be  taken.  Using  Inequality  7  is  the  more 
systematic  method.  The  educated  guess  method,  however,  illus¬ 
trates  the  effect  of  4  on  adaptive  convergence  time. 

An  educated  guess  of  4  =  0.1  was  made.  With  4  =  0.1,  the 
canceler  output  for  N  -  16,  32  64  and  128  converged.  The  output 

for  N  =  256  did  not  converge.  When  4  was  set  equal  to  0.04  the 
canceler  output  for  N  =  256  converged  quite  rapidly. 

The  simulation  outputs  for  N  =  16,  32,  64  and  128  and  4  = 

0.1  are  shown  in  Figures  6  (a,b,c),  7  (a,b,c),  8  (a,b,c),  and 
9  (a,b,c),  respectively.  The  outputs  for  N  =  256  and  4  =  0.04 
■'•■re  shown  in  Figure  10  (a,b,c).  The  "a"  figure  of  each  set  shows 
the  input  parameters  and  the  canceler  power  output  for  each 
iteration.  The  iterations  are  continued  until  the  user  specified 
reduction  in  interfering  signal  strength  (30  dB)  has  been 
achieved.  Notice  that  as  the  given  value  of  4  (=0.1)  gets  closer 
a  .he  optimum  value  of  4  (by  increasing  N  in  Inequality  7)  the 
.amber  of  iterations  necessary  to  achieve  the  user-specified 
reduction  in  interfering  power  decreases.  Once  the  time  neces- 
■ary  to  complete  a  single  iteration  is  known,  or  assumed,  the 
^ctal  time  necessary  to  reduce  the  interfering  signal  the  re- 
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quired  number  of  dB,  the  "adaptive  convergence  time",  can  be 

calculated . 

The  "a"  figure?  give  the  canceler  output  power  versus 
iteration  number,  but  they  tell  nothing  about  output  power  versus 
frequency.  To  work  properly  the  canceler  must  attenuate  the 
interfere!',  but  not  the  intended  signal.  The  PTF  must,  there¬ 
fore,  pass  the  interferer  with  a  zero  phase  shift  and  attenuate 
the  intended  signal.  PTF  amplitude  gain  and  phase  shift  are 
given  in  the  "b"  and  "c"  figures,  respectively,  of  each  set. 

The  "b"  figures  show  that  as  N  increases,  the  width  of  the 
central  lobe  of  the  PTF  amplitude  gain  curve  decreases.  PTF 
frequency  resolution  has  been  previously  defined  as  the  differ¬ 
ence  in  frequency  between  the  center  of  the  PTF  amplitude  gain 
curve  and  its  first  zero.  Resolution  is  given  by  Equation  6. 

1  1 

Frequency  Resolution  =  -  = - (6) 

NA  (total  delay) 

where : 

N  =  number  of  taps 
A  =  intertap  delay 

As  N  increases,  the  frequency  resolution  of  the  PTF  gets 
better,  i.e.,  the  minimum  frequency  separation  needed  for  the  PTF 
to  pass  a  strong  interferer  and  reject  a  weak  intended  signal, 
decreases.  Figures  6b,  7b,  8b,  9b,  and  10b  clearly  show  this 
relationship  between  N,  the  number  of  PTF  taps,  and  frequency 
resolution.  For  each  curve,  the  frequency  separation  between  the 
curve's  center  frequency  and  first  zero  is  given  by  Equation  6. 
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(text  continues  on  page  50) 


run 

Enter  the  number  of  taps  in  the  transversal  filter  16 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.1 

Enter  intended  signal  frequency  (in  MHz)  ISO 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  1S1 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz>  150 
Enter  the  highest  expected  frequency ( in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 
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Figure  6a.  Iterated  Canceler  Output  Power  for  N  =  16  taps 
MU  =0.1. 
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PTF  Amplitude  Gain  us  Frequency 

weights  calculated  uia  Midrows  algorithm 


re  6b.  PTF  Amplitude  Gain  vs  Frequency  for  N  =  16  taps 


PTF  Phase  Shift  ms  Frequency 

weights  calculated  uia  Widrows  algorithm 


Figure  6c.  PTF  Phase  Shift  vs  Frequency  for  N  =  16  taps 
MU  =  0.1. 


r  vn 

Enter  the  number  of  taps  in  the  transversal  filter  32 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.3444 

Enter  MU  the  convergence  parameter  .1 

Enter  intended  signal  frequency  (in  MHz)  1S0 

Enter  intended  signal  strength  (in  dBm)  -95 

Enter  interfering  frequency  (in  MHz)  151 

Enter  interfering  signal  strength  (m  dBmJ  20 

Enter  desired  reduction  (in  d6)  of  interfering  signal  strength  30 
Enter  3IMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequencyfin  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  dbJ  0 


Iteration  Output  Power  (dBm) 
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Figure  7a.  Iterated  Canceler  Output  Power  for  N  s  32  Taps, 
MU  =  0.1. 
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RTF  Amplitude  Gain  us  Frequency 

weights  calculated  uia  Midrows  algorithm 


Figure  7b.  PTF  Amplitude  Gain  vs  Frequency  for  N  =  32  taps 
MU  =  0.1. 


run 

Enter  the  number  of  taps  in  the  transversal  filter  64 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  .1 

Enter  intended  signal  frequency  (in  MHz)  1S0 

Enter  intended  signal  strength  (in  dBm)  -93 

Enter  interfering  frequency  (in  MHz)  131 

Enter  interfering  signal  strength  (in  dBm;  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  3IMTIME  (a  dimens lonless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz>  150 
Enter  the  highest  expected  frequency(in  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  db)  0 
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Figure  8a.  Iterated  Canceler  Output  Power  for  N  s  64  Taps, 
MU  =  0.1. 
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PTF  Amplitude  Gain  os  Frequency 

weights  calculated  via  Widrows  algorithm 
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Frequency  (Mhz) 

Figure  8b.  PTF  Amplitude  Gain  vs  Frequency  for  N  =  64  taps 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 
nter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 
nter  MU  the  convergence  parameter  .1 
Enter  intended  signal  frequency  (in  MHz)  180 
Enter  intended  signal  strength  (in  dBm)  -93 
Enter  interfering  frequency  (in  MHz)  131 
Enter  interfering  signal  strength  (in  dBm)  20 

Enter  desired  •'eduction  (in  dB)  of  interfering  signal  strength  30 
Enter  3IMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency(in  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  db)  0 
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Figure  9a.  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 
MU  =  0.1. 
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F’TF  Amplitude  Gain  os  Frequency 

weights  calculated  oia  Widrows  algorithm 
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Figure  9b.  PTF  Amplitude  Gain  vs  Frequency  for  N  =  128  taps 
MU  =0.1. 


PTF  Phase  Shift  us  Frequency 


0.  £  fl)  i.H  li  >A  JZ  <*■  •+»  C  73  <Xi  !Ti  S_  <1)  1»  ifi 


46 


Figure  9c.  PTF  Phase  Shift  vs  Frequency  for  N  =  128  taps 


r  un 

Enter  the  number  of  taps  1  n  the  transversal  filter  £56 

Enter  the  del an  tine  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.04 

Enter  intended  signal  frequency  (in  MHz)  ISO 

En  ter  i n  tended  s l gna 1  s tr eng  th  ( i n  dBm )  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz>  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTE  voltage  loss  per  tap  (in  db)  0 

Iteration  Output  Power  (dBm) 

1  1.70E+-1  OPTIMUM  MU=0. 0390625 

2  -1.69E+1 


Figure  10a.  Iterated  Canceler  Output  Power  for  N  =  256  taps 
MU  =  0.04. 
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FTP  Amplitude  Gain  us  Frequency 

weights  calculated  uia  Midrows  algorithm 


Figure  10b.  PTF  Amplitude  Gain  vs  Frequency  for  N  =  256  taps 
MU  =  0.04. 


PTF  Phase  Shift  ms  Frequency 


Frequency  (Mhz) 

Figure  10c.  PTF  Phase  Shift  vs  Frequency  for  N  =  256  taps 
MU  =  0.04. 


Therefore,  in  order  to  choose  an  appropriate  value  of  N, 
either  use  the  required  frequency  resolution  to  solve  for  N  (via 
Equation  6)  or  use  the  PTF  amplitude  gain  vs.  frequency  curve 
generated  by  the  simulation  to  give  a  more  detailed  look  at  the 
frequency  transfer  function  (as  a  function  of  N) . 

For  N  =  128,  144  and  256,  the  respective  frequency  resolu¬ 
tions  (via  Equation  7)  are  1.125  MHz,  1  MHz  and  0.506  MHz. 

Setting  N  =  144  or  256  both  achieves  or  betters  the  required 
resolution  of  1  MHz.  But  they  require  more  complicated  and 
expensive  PTFs  and  associated  circuitry  than  for  N  =  128,  simply 
because  they  have  more  taps.  N  was  chosen  equal  to  128  because: 
(1)  it  almost  achieves  the  "required  resolution"  of  1  MHz,  (2)  a 
128  tap  PTF  will  be  somewhat  simpler  than  a  144  or  256  tap  filter 
and,  (3)  we  are  interested  in  modeling  a  Texas  Instruments  PTF 
being  developed  under  Contract  DAAL01-88-C-0831  that  has  128 
taps . 

Figures  6c,  7c,  8c,  9c,  and  10c  show  that  the  PTF  phase 
shift  versus  frequency  curve  fluctuates  faster  as  N  increases. 

The  PTF  phase  shift  must  behave  in  this  manner  in  order  for  the 
canceler  to  be  able  to  separate  signals  progressively  closer  in 
frequency.  An  adaptive  noise  canceler  works  by  having  the  PTF 
pass  the  strong  interfering  frequency  with  a  zero  phase  shift  (so 
that  the  PTF  filtered  interferer  subtracts  in  phase  at  the  summer 
of  the  canceler  from  the  non-filtered  interferer) .  In  addition, 
it  attenuates  the  weak  intended  signal  with  a  non-zero  phase 
shift  (so  that  the  PTF  filtered  intended  signal  subtracts  out  of 
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phase  at  the  summer  from  the  non-filtered  intended  signal,  thus 
minimizing  the  effect  of  this  subtraction  on  the  unfiltered  input 
signal).  The  closer  in  frequency  that  the  interfering  and 
intended  signals  are,  the  faster  the  PTF  phase  shift  will  have  to 
change  in  order  for  there  to  be  a  significant  difference  in  the 
phase  shift  at  the  two  frequencies.  N  has  to  be  large  enough  so 
that  the  phase  shift  varies  fast  enough  to  insure  signal  separa¬ 
tion  . 

Since  N  is  now  fixed  equal  to  128,  the  only  design  parameter 
left  to  be  determined  is  the  convergence  parameter  /x. 

With  an  educated  guess  of  /n  =  0.1,  four  iterations  were 
necessary  to  reduce  the  interference  by  30  dB.  A  value  of  /x  that 
reduces  the  number  of  iterations  below  four  is  needed  in  order  to 
minimize  the  adaptive  convergence  time.  Two  iterations  is  the 
minimum  number  of  iterations  necessary  for  significant  interfer¬ 
ence  reduction.  Significant  interference  reduction  cannot  be 
achieved  in  one  iteration,  since  in  the  first  iteration  the 
simulation  sets  all  tap  weights  equal  to  zero.  As  a  result, 
there  is  only  a  3  dB  reduction  in  power  due  to  the  action  of  the 
power  splitter  in  the  adaptive  noise  canceler.  In  other  words, 
in  Figure  3,  assuming  all  the  tap  weights  are  equal  to  zero 
inplies  there  is  no  output  from  the  PTF.  In  effect,  half  the 
input  power  is  being  lost  and  there  is  a  3  dB  reduction  in 
canceler  output  power  relative  to  the  input  power. 

Tabic  1  summarizes  the  results  of  a  number  of  simulation 
runs  with  different  values  of  /i .  The  number  of  iterations 
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iic.ess.iry  co  reduce  the  interference  by  30  dB  is  a  very  sensitive 
function  of  /i .  in  the  approximate  region  of  /x  -  0.075  to  /x  = 
0.03,  30  dB  interference  reduction  is  achieved  in  2  iterations. 
This  region  corresponds  to  the  upper  bound  of  Inequality  7. 

1  1  1 

0  v  jx  < - =  -  = - 

(N  +  1)  (input  power)  (128  +  1)  20  dBm  129  (0.1  watt) 

So  the  upper  bound  of  Inequality  7  can  be  used  to  derive  an 
optimal  value  (or  range  of  values)  for  /x. 

From  Table  1,  the  optimal  range  for  /x  is  from  0.075  to 
0.08.  The  convergence  parameter  will  "arbitrarily"  be  assumed  to 
be  0.08.  The  design  of  the  adaptive  filter  is  now  complete. 

An  adaptive  noise  canceler  with  design  parameters  assumed 
here  (number  of  PTF  taps  =  128,  intertap  delay  =  6.944  nanosec¬ 

onds,  and  n  =  0.08)  will  reduce  interfering  or  strong  signals  by 
30  dB  in  two  iterations  while  having  a  minimal  affect  on  weak 
intended  signals  1.125  MHz  away  from  the  interferer. 

Figures  20a  and  20b  give  the  PTF  amplitude  gain  vs.  frequ¬ 
ency  and  phase  shift  vs.  frequency  respectively  for  N  =  128  and 
M  =  0.08.  Notice  that  these  figures  are  identical  to  figures  9b 
and  9c  for  N  =  128  and  /x  =  0.1.  This  is  so  because  ju  does  not 
affect  the  optimal  PTF  frequency  transfer  function.  It  only 
affects  the  number  of  iterations  necessary  to  reach  it. 

It  should  be  re-emphasized  that  Inequality  7  implies  that  if 
;x  is  chosen  within  the  range  of  the  Inequality  then  the  expected 
weight  vector  will  converge  to  the  optimal  or  Wiener  weight 
vector.  It  should  not  necessarily  converge  in  a  minimum  number 
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Table  1.  Convergence  Parameter  vs  Number  of  Iterations 

Necessary  for  a  30  db  Reduction  in  Interference 
Power 


Number  of  Iterations  Necessary  for  30  db 
Reduction  in  Interference  power 

Figure 

0 . 2 

No  convergence  -  simulation  output  "blows  up" 

11 

0. 16 

No  convergence  -  output  oscillates 

12 

0. 15 

36 

13 

0.1 

4 

9a 

0.0 

3 

14 

CO 

o 

o 

2 

15 

1 

0 . 0775 

2 

16 

0.075 

2 

17 

0.07 

3 

18 

!  0.06 

4 

19 

5  3 


of  iterations.  Yet,  this  is  exactly  what  the  simulation  runs 
indicate.  If  jx  =  1/ ( (N+l) ( interfering  power))  then  the  adaptive 
process  converges  in  two  iterations.  The  reason  for  this  is  not 
absolutely  clear.  What  is  probably  happening  is  that  when 
a  --  1/ (  (N+l )( interfering  power)),  the  adaptive  process  is  criti¬ 
cally  damped.13  Critical  damping  gives  convergence  in  one 
iteration  (two  iterations  in  our  simulation  because  the  initial 
weight  vector  is  assumed  equal  to  zero  and  we  are  counting  it  as 
an  iteration).  In  addition,  if  /x  is  greater  than  twice  the  /x 
value  that  gives  critical  damping,  then  the  adaptive  process  does 
not  converge.  This  is  exactly  what  Table  1  shows. 
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(text  continues  on  page  69) 


Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.2 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

En ter  i n  ter  t  er i ng  f r equenc y  ( i n  MHz )  181 

Enter  interfering  signal  strength  (in  dBm)  20 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequencytin  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  db)  0 

I  ter  a  t i on  Qu  tpu t  Power  ( dBm ) 


1 

1  .70E+1 

2 

2 .08E+1 

~J 

9 .22E+0 

4 

2.02E+1 

tr 

_i 

3.12E+1 

K 

4.22E+1 

7 

5.33E+1 

8 

6.43E+1 

9 

7.53E+1 

10 

8.63E+1 

1  1 

9.73E+1 

12 

1 .08E+2 

Figure  11.  Iterated  Canceler  Output  Power  for  N  =  128  taps 

MU  =  0.2. 
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1 .19E+3 
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38 
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31 
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33 

3.40E+3 
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3 . 63E+3 
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37 

3.84E+3 
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Figure  11.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.2.  (continued) 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.16 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SI NT I ME  (a  dimensionless  parameter  between  8  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency ( in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  F’TF  voltage  loss  per  tap  (in  db)  0 


Iteration  Output  Power  (dBm) 


1 

1  .78E+1 

£ 

1  .73E+1 

1  .6£E+1 

4 

1  . 49E+1 

cr 

1  .  5£E+1 

K 

1  .56E+1 

r 

1  .59E+1 

o 

1 .63E+1 

1  .66E+1 

10 

1  .70E+1 

Figure  12.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.16. 
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1 .6£E+1 

13 

1 .66E+1 

14 
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1 .67E+1 
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1 .65E+1 
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1 .66E+1 

£6 

1 .6£E+1 
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U.  1 

1 .66E+1 

£8 

1 .69E+1 

£9 

1 .67E+1 
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1 .65E+1 
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1 .63E+1 
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1 .7£E+1 
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1 .64E+1 
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1 .59E+1 
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1 .6£E+1 

37 
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1 .69E+1 

39 

1  .67E+1 

Figure  12.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.16.  (continued) 
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run 

Enter  the  number  of  taps  in  the  transuersal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  .15 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  20 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency(in  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap. (in  db)  0 


Iteration  Output  Power  (dBm) 


1 

1 .70E+1 

2 

1  .62E+-1 

r> 

1  .54EM 

4 

1 .46E+1 

C* 

1 . 39E+1 

6 

1  .31E+1 

7 

i 

1 .23E+1 

3 

1  . 15E+1 

o 

1 .07E+1 

10 

9.96E+0 

Figure  13.  Iterated  Canceler  Output  Power  for  N  =  128  Tape, 
MU  =  0.15. 
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36 
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Figure  13.  Iterated  Canceler  Output  Power  for  N  =  128  Taps, 
MU  =  0.15.  (Continued) 
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run 

Enter  the  number  of  taps  in  the  transuersal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  .09 

Enter  intended  signal  frequency  (in  MHz)  1S0 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  20 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  3IMTIME  (a  dimensionless  parameter  between  0  and  4)  i 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency(in  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  db)  0 


Iteration  Output  Power  (dBm) 

1  1.70E+1 

2  4.23E-1 

3  -1.60E+1 


Figure  14.  Iterated  Canceler  Output  Power  for  N  a  128  Taps, 
MU  =  0.09. 
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r  un 

Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds.!  6.9444 

Enter  MIJ  the  convergence  parameter  .08 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Erter  SIMTIME  (a  dimensionless  parameter-  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 


Iteration  Output  Power  (dBm) 


1  . 70E+1 
-1 .65E+1 


Figure  15.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.08. 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  del  an  time  between  taps  (in  nanoseconds.1  6.9444 

Enter  MU  the  convergence  parameter  0.0775 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  PTF  voltage  loss  per  tap  (in  db)  0 

Iteration  Output  Power  (dBm) 

1  1.70E+1 

£  -2.16E+1 


Figure  16.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.775. 
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r  u  r  i 

Enter  the  number  of  taps  in  the  transversal  filter  1£8 

Enter  the  delay  time  between  tans  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  ©.075 

Enter  intended  signal  frequency  (in  MHz)  1S6 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 


I  ter a t i on 


Output  Power  (dBm) 


1  1.70E+1 

£  -1 .83E+1 


Figure  17.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.075. 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 
Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 
Enter  MU  the  convergence  parameter  9.0? 

Enter  intended  signal  frequency  (in  MHz)  180 
Enter  intended  signal  strength  (in  dBm)  -98 
Enter  interfering  frequency  (in  MHz)  181 
Enter  interfering  signal  strength  (in  dBm)  £0 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 

Iteration  Output  Power  (dBm) 

1  1.70E+1 

£  -2.43E+0 

3  -2.14E+1 


Figure  18.  Iterated  Canceler  Output  Power  for  N  =  128  taps 
MU  =  0.07. 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.66 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  20 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  210 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 

Iteration  Output  Power  (dBm) 

1  1 .70E+1 

2  4.39E+0 

3  -3.19E+0 

4  -2.04E+1 


Figure  19.  iterated  Canceier  Output  Power  for  N  =  128  taps 
MU  =  0.06. 
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PTF  Amplitude  Gain  us  Frequency 

weights  calculated  via  Widrows  algorithm 
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Amplitude  Gain  vs  Frequency  for  N  =  128 


PTF  Phase  Shift  vs  Frequency 

weights  calculated  via  Uidrows  algorithm 


Frequency  (Mhz) 

PTF  Phasa  Shift  vs  Frequency  for  N  =  128 
Taps,  MU  =  0.08. 


The  design  example  has  assumed  that  the  input  interfering 
power  has  been  limited  to  +20  dBm  by  a  frequency  selective 
limiter.  Since,  in  our  example,  the  input  interfering  power  is 
always  greater  than  or  equal  to  +20  dBm,  the  output  of  the 
frequency  selective  limiter  is  a  constant  +20  dBm.  This  is  the 
constant  power  that  the  adaptive  noise  canceler  "sees." 

Adaptive  convergence  time,  i.e.,  the  number  of  iterations 
necessary  to  reduce  the  interfering  power  by  a  user  given  amount, 
is  a  function  of  the  interfering  power  level.  "Large"  interfer- 
ers  can  be  reduced  a  fixed  amount  (e.g.,  30  dB)  faster  (i.e.,  a 
smaller  number  of  iterations)  than  "small"  interferers. 

This  is  illustrated  in  Figure  21.  Figure  21  shows  a  simula¬ 
tion  run  in  which  the  design  parameters  of  Figure  15  (/x  =  0.08,  N 
=  128,  intertap  delay  =  6.9444  nanoseconds)  are  used  on  an 
interfering  signal  that  is  10  dB  less  than  the  interfering  signal 
power  used  in  Figure  15  (10  dBm  vs.  20  dBm). 

The  output  in  Figure  15  dropped  over  30  dB  in  2  iterations. 
The  output  in  Figure  21  took  30  iterations  to  drop  30  dB.  This 
behavior  can  be  explained  using  Equation  1,  Widrow's  LMS  algo- 
r ithm: 

— ¥  -* 

wk+l  =  wk  +  2^ekxk  (*) 

»here : 

Wj^+  ^  =  tap  weight  vector  at  the  k+lth  iteration. 

Wj,  -  tap  weight  vector  at  the  kth  iteration. 

v  =  convergence  parameter. 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  128 

Enter  the  delay  time  between  taps  (in  nanoseconds.1  6.9444 

Enter  MU  the  convergence  parameter  .08 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  10 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  SIMTIME  i.  a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  frequency (in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 

Iteration  Uutput  F'ower  (dBm) 


1 

6 . 99E+0 

c 

6 .05E+0 

■g 

5.12E+0 

4 

4.13E+0 

cr 

3.25E+0 

if, 

£ . 31E+0 

i 

1 . 33E+0 

:f; 

4.43E-1 

q 

-4.92E-1 

10 

-1 . 43E+0 

Figure  21 .  Iterated  Canceler  Output  Power  for  N  =  128  taps 

MU  =  0.08,  and  Interfering  Signal  Strength  =  10  dBm. 
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-1 .8£E+1 
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-1  .9£E+1 

30 

-£ .01E+1 

Figure  21 .  Iterated  Canceler  Output  Power  for  N  =  128  taps 

MU  =  0.08,  and  Interfering  Signal  Strength  =  10  dBm. 
(continued) 
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dune -.'I  :T  err.-'  utput  at  the  kth  interation 


X ^  PTF  tap  amplitude  vector  at  the  kth  iteration. 


When  the  input  interfering  power  is  reduced  both  and  Xk 
will  be  reduced.  If  the  input  interfering  power  is  reduced  by  a 
factor  of  10,  then,  since  both  ek  and  Xk  are  amplitudes,  they 
will  be  reduced  by  JlO  and  their  product  will  be  reduced  by  a 
factor  of  10.  As  a  result,  the  amount  by  which  the  weights  are 
incremented  at  each  iteration,  2  g  ek  X^,  is  decreased  by  a 
factor  of  10.  It  now  takes  the  adaptive  process  much  longer  to 
achieve  the  required  power  reduction. 

Equation  1  also  suggests  a  way  of  remedying  the  problem. 

—v 

Simply  increase  g  in  the  same  proportion  that  X ^  was  de¬ 
creased,  i.e.,  by  a  factor  of  10.  This  was  done  in  a  simulation 
run  illustrated  in  Figure  22  in  which  g  was  increased  from  g  = 
0.08  (in  Figure  21)  to  g  =  0.8.  The  output  drops  over  30  dB  in 
just  two  iterations.  The  "problem'*  has  been  fixed. 

Since  input  interfering  power  affects  the  adaptive  conver¬ 
gence  time,  then,  in  any  simulation  of  a  frequency  hopping  cosite 
problem,  either  the  adaptive  noise  canceler  must  "see"  a  suffi¬ 
ciently  high  interference  power  to  minimize  the  adaptive  conver¬ 
gence  time  (for  a  given  g)  or  u  must  be  recalculated  via  Equation 


7 


1 

u  ~  - 

(N  +  1)  (interfering  power) 


(7) 


every  time  the  input  power  is  measured. 
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run 

Enter  the  number  of  taps  in  the  transversal  filter  1£S 

Enter  the  delay  time  between  taps  (in  nanoseconds)  6.9444 

Enter  MU  the  convergence  parameter  0.8 

Enter  intended  signal  frequency  (in  MHz)  180 

Enter  intended  signal  strength  (in  dBm)  -98 

Enter  interfering  frequency  (in  MHz)  181 

Enter  interfering  signal  strength  (in  dBm)  10 

Enter  desired  reduction  (in  dB)  of  interfering  signal  strength  30 
Enter  3IMTIME  (a  dimensionless  parameter  between  0  and  4)  1 
Enter  the  lowest  expected  signal  frequency  (in  MHz)  150 
Enter  the  highest  expected  f'requency(in  MHz)  £10 

Enter  the  frequency  increment  (in  KHz)  to  be  used  for  plotting  100 
Enter  the  RTF  voltage  loss  per  tap  (in  db)  0 

Iteration  Output  Power  (dBm) 

1  6.99E+9 

£  -£.65E+1 


Figure  22.  Iterated  Canceler  Output  Power  for  N  =  128  taps 

MU  =  0.8,  and  Interfering  Signal  Strength  =  10  dBm. 
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Up  to  this  point  the  PTF  was  assumed  to  have  no  tap  loss,  as 
discussed  in  the  "Simulation  Input  Parameters"  section.  Once  the 
tap  loss  is  known,  then  the  optimal  /i  determined  for  zero  tap 
loss  should  be  multiplied  by  (1/tap  loss)  to  get  the  actual  value 
of  n  to  be  used  when  building  a  real  adaptive  noise  canceler. 
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CONCLUSION 


The  adaptive  noise  canceler  simulation  is  able  to  handle  the 
typical  parameters  encountered  when  two  frequency  hopping  radios 
are  co-located.  It  can  identify  design  parameters  that  cause 
the  canceler  to  reduce  a  +20  dBm  interfering  signal  by  30  dB  in 
one  or  two  iterations,  depending  on  the  definition  of  an  iteration. 
The  simulation  is  suitable  for  use  as  an  adaptive  noise  canceler 
design  tool  to  evaluate  the  effect  of  design  parameter  changes  on 
canceler  performance  by  determining  the  adaptive  convergence  time 
and  PTF  frequency  transfer  function. 
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